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We have searched for diphoton events {jjI^t) with large missing trans- 
verse momentum, ^J^t events {■jI^t+ ^ 2 jets) with two or more 
jets, and diphoton events (77) with high transverse energies in pp col- 
hsions at ^/s = 1.8 TeV using approximately 100 pb~^ of data col- 
lected with the D0 detector at the Fermilab Tevatron in 1992-1996. 
No excess of events beyond the expected backgrounds is observed. The 
null results are interpreted in supersymmetric models with a dominant 
X2 ~^ 7X1 decay and in terms of Dirac pointlike monopole production. 
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in Particle Physics, La Thuile, Italy, March 1-7, 1998. 
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I. INTRODUCTION 



The D0 detector collected a data sample corresponding to an integrated luminosity of 
approximately 100 pb^^ during the 1992-1996 Tevatron Run at ^/s = 1.8 GeV. The detector 
consists of three major components: a non-magnetic tracking system including a transition 
radiation detector in the central region, a liquid-argon calorimeter, and a toroidal magnetic 
muon spectrometer. A detailed description of the D0 detector can be found in Ref. 
The Tevatron pp collider is ideal for studying high p^ phenomena. It is a natural place to 
search for new physics beyond the Standard Model. In this talk, the results of the searches 
for supersymmetry and Dirac monopoles are presented. 

II. SEARCH FOR SUPERSYMMETRY USING PHOTONS 

Supersymmetry P| is a generalization of space-time symmetry. It introduces for every 
Standard Model particle a supersymmetric partner that differs in spin by 1/2. 7^-parity con- 
servation requires that supersymmetric particles be produced in pairs and that the lightest 
supersymmetric particle (LSP) be stable. 

The minimal supersymmetric standard model (MSSM) is the simplest supersymmetric 
model. In the MSSM, the Gaugino-Higgsino sector (excluding gluinos) is described by four 
parameters: Mi, M2, fi, and tan/5, where Mi and M2 are the U{1) and SU{2) gaugino 
mass parameters, n is the Higgsino mass parameter, and tan f3 is the ratio of the vacuum 
expectation values of the two Higgs doublets. The Gaugino-Higgsino mixing gives in four 
neutral mass eigenstates (neutralinos Xi*; i = 1;2,3,4) and two charged mass eigenstates 
(charginos x^, j = 1, 2) whose masses and couplings are fixed by the above four parameters. 
Supersymmetric models |^,|] predicting photon production have been proposed as possible 
explanations of a recent event [^] reported by the CDF collaboration. Within the MSSM, it 
has been shown that the radiative decay of X2 IXi dominates in some regions of parameter 
space 0. Assuming x1 is the LSP, the production of X2, either directly or indirectly from 
decays of other supersymmetry particles, will result in events with two high transverse energy 
[Et) photons and large missing transverse momentum (I^t) and/or in events (7^7^ + ^ jets) 
with one high Et photon, multile jets and large I^t- Moreover, the '~f'~fI^T events are also 
expected in supersymmetric models with a light gravitino (G) being the LSP. In this case, 
Xi (assumed to be the lightest superpartner of a standard model particle) is unstable and 
decays into a photon plus a gravitino (x? iG). 

We present searches for new physics in the channels pp — >■ + X and pp 7-^t+ ^ 

2 jets -|- X at the Fermilab Tevatron. The 7^r events with fewer than two jets are not 
considered here, due to large backgrounds from H^-|-jets production. We interpret our results 
of the 77-^T analysis in terms of ee, z>i>, and X2X2 production and of the '~iIfT+ > 2 jets 
analysis in terms of squark (g) and gluino [g) production in supersymmetric models with a 
dominant x\ ~^ IXi decay. The interpretation of the results of the jjI^t analysis in models 
with a light G can be found in Ref. 0. 

The trigger used requires one electromagnetic (EM) cluster with transverse energy Ex > 
15 GeV, one jet with Et > 10 GeV, and ^t> 14 GeV (^t> 10 GeV for about 10% of 
the data taken early in the Tevatron run). The jets in the trigger include non-leading EM 
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clusters. Photons are identified via a two-step process: the selection of isolated EM energy 
clusters and the rejection of electrons. The EM clusters are selected from calorimeter energy 
clusters by requiring (i) at least 95% of the energy to be deposited in the EM section of the 
calorimeter, (ii) the transverse and longitudinal shower profiles to be consistent with those 
expected for an EM shower, and (iii) the energy in an annular isolation cone with radius 0.2 
to 0.4 around the cluster in r/ — space to be less than 10% of the cluster energy, where r) 
and (f) are the pseudorapidity and azimuthal angle. Electrons are removed by rejecting EM 
clusters which have either a reconstructed track or a large number of tracking chamber hits 
in a road between the calorimeter cluster and the event vertex. I^t is determined from the 
energy deposition in the calorimeter for |?7| < 4.5. 

A. Search for Events 

To be selected as candidates, events are first required to have two identified pho- 

tons, one with E^p > 20 GeV and the other with Ej? > 12 GeV, each with pseudorapidity 
I'q'^l < 1.1 or 1.5 < \ri'^\ < 2.0, the regions with good photon identification. We denote the 28 
events passing these photon requirements as the 77 sample. We then require I^t> 25 GeV 
with at least one reconstructed vertex in the event to ensure good measurement of I^t- No 
requirement on jets is made. Two events satisfy all requirements. The data used in this 
analysis correspond to an integrated luminosity of 106.3 ± 5.6 pb~^. 

The principal backgrounds are multijet, direct photon, W + •y, W + jets, Z ee, and 
Z ^ TT ^ ee events from Standard Model processes with misidentified photons and/or 
mismeasured ^t- The background due to I^t mismeasurement is estimated using events 
with two EM-like clusters which satisfy looser EM cluster requirements than those discussed 
above, and for which at least one of the two fails the EM shower profile consistency require- 
ment (ii) above. In addition, these events must pass the photon kinematic requirements. 
By normalizing the number of events with < 20 GeV in the QCD sample to that in the 
77 sample, we obtain a background of 2.1 ± 0.9 events due to I^t mismeasurement for 1^^ 
> 25 GeV. 

Other backgrounds are due to events with genuine I^t such as those from W+^'y' (where 
'7' can be a real or a fake photon), Z ^ tt —>■ ee, and ti — >■ ee + jets production. These 
events would fake yylpT events if the electrons were misidentified as photons. We estimate 
their contribution using a sample of e'7' events passing the kinematic requirements, including 
that on Ifx- Taking into account the probability (0.0045 ± 0.0008, determined from Z — > ee 
data) that an electron is misidentified as a photon, we estimate a background of 0.2 ± 0.1 
events. Adding the two background contributions together yields 2.3 ±0.9 events. No excess 
of events is observed. 

The null results are used to set upper cross section limits on the pair production of 
scalar electron (e), scalar neutrino (i>), and X2 with Br(e — > ex^) = 100%, Br(z/ — ^ z/^^) = 
100%, and Br(x2 7Xi) = 100%. pp ee, X2X2 events are generated using the 
ISAJET program and are processed through the detector and trigger simulation and the 
reconstruction program. The e mass is set to 100 GeV/c^ while the z/, Xi masses 

are varied between 70-100 GeV/c^, 50-90 GeV/c^, and 30-80 GeV/c^ respectively. The 
efficiency is found to be insensitive to the e and 9 masses and is largely determined by 
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(a) The total efficiency and (b) the 95% C.L. cross section hmit as a function of m^o—m^o 



FIG. 1. 

for pp — > 
systematic error is estimated to be 



ee, i>D, X2X2 processes. 



The error shown on the efficiency is statistical only and the 



the mass difference m^o — m^o as shown in Fig. 0(a). For a given value of m^o — m^o, 
the efficiencies for the three processes are the same within the errors, independent of the 
values of mg, m^, and m^o. The 95% confidence level (C.L.) upper limits on the cross 
section as a function of m^o — m^o is shown in Fig. ^(b). The limit is about 300 fb for large 
values of m^o — m^o. These results improve over those we published earlier [Q, mainly due 
to the improved photon identification. However, we note that the expected cross sections 
for pp ee, uu, X2X2 processes are small, well below the experimental upper limits. For 
example, the theoretical cross section for pp eiei + e.R(iR production is about 30 fb for 
me = 100 GeV/c2. 



B. Search for j^t + n Jets Events 

To be selected as 'yI^T+ > 2 jets candidates, events are first required to have at least 
one identified photon with Ej. > 20 GeV and pseudorapidity \ri'^\ < 1.1 or 1.5 < \ri'^\ < 2.0, 
and two or more jets with Ej^ > 20 GeV and \r]^ < 2.0. We denote events passing these 
requirements as the 7+ > 2 jets sample. The I^t distribution of the 7+ > 2 jets events is 
shown in Fig. H(a). After requiring I^t > 25 GeV, 378 events remain in the sample. The 
data used in this analysis correspond to an integrated luminosity of 99.4 ± 5.4 pb~^. 

The principal backgrounds are multijet, direct photon, VT+jets, and Z+jets events from 
Standard Model processes with jets or electrons misidentified as photons and/or mismea- 
sured I^t- Following the same procedure as in the 77^t analysis, the background due to 
I^T mismeasurement is estimated to be 370 ± 36 events. W+ > 2 jets events with W ^ 
would fake 'yI^T+ > 2 jets events if the electrons were misidentified as photons. This con- 
tribution is estimated to be 4 ± 1 using a sample of e^r+ > 2 jets events passing all the 
kinematic requirements with the electron satisfying those imposed on the photon. Another 
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background due to iu) + jets and Z{-^ z/z/) + jets production is found to be negligible. 

The total background (374 ± 36) agrees well with the number of observed events. The Ht 
(defined as the scalar sum of the Et of all jets with Et > 20 GeV and |?7| < 2.0) distribution 
is shown in Fig. ||(b) for both ■jI^t+ > 2 jets and background samples. The background 
distribution reproduces the observed 'jI^t+ > 2 jets distribution well. 
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FIG. 2. (a) The I^t distribution (solid circles) of the events with one photon and two or more jets. 
The expected I^t distribution from the backgrounds is shown as a solid histogram. The number 
of events with 1^t< 20 GeV in the background is normalized to that in the 7+ > 2 jets sample, 
(b) The Ht distributions of the ^Ij^T+ > 2 jets and background events. Also shown (dashed and 
dotted histograms) are the expected distributions from supersymmetry for two different values 
squark/gluino masses assuming an equal squark and gluino mass. 



With the basic selection discussed above, the 7-^t+ > 2 jets sample is dominated by 
the backgrounds. To optimize the selection criteria for supersymmetry, we simulate squark 
and gluino production (either in pair or in association with charginos/neutralinos) using 
the Pythia program [0]. The MSSM parameter values are set to Mi = M2 = 60.0 GeV, 
tan/5 = 2.0, and /i = —40.0 GeV. This set of parameter values gives m^o = 33.5 GeV/c^, 

m^o = 60.0 GeV/c^, and Br(x2 7Xi) = 100%. Sleptons (£) and stop (i) are assumed to be 
heavy. Events with X2 f he final state are selected and are run through a detector simulation 
program, a trigger simulator, and the same trigger requirements, reconstruction, and analysis 
as the data. The I^t and Ht distributions for mq(= rrig) = 200,300 GeV/c^ events after 
the basic selection are shown in Figs. 0. The distributions expected from supersymmetry 
are considerably harder than those of the background events. To increase sensitivity to 
supersymmetry, we introduce an Ht cut and maximize the e/ab ratio by varying ^t and 
Ht requirements. Here e is the efficiency for mq{= rrig) = 300 GeV/c^ MC events and ab is 
the error on the estimated number of background events. The optimized cuts are found to 
be IpT > 45 GeV and Ht > 220 GeV. With these additional cuts, five 7-^t+ > 2 jets events 
are observed while 8 ± 6 events are expected from the background processes. 

The efficiency for supersymmetry signal varies from a few percent for low mass q/g to 
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approximately 25% for high mass q/g. MC studies show that the efficiency varies by 4% for 
different choices of Mi, M2, tan/3, and /i which are consistent with Br(x2 7Xi) = 100% 
and m^o — m^o > 20 GeV/c^ suggested in Ref. 0]. Experimentally, the mass requirement 
is needed to ensure that photons from X2 decays are reasonably energetic and are detected 
with good efficiency. The variation is assigned as a systematic error in the efficiency. The 
total fractional systematic error on the efficiency is 9%. 

With five events observed and 8 ± 6 events expected from backgrounds, we observe 
no excess of events. We compute 95% C.L. upper limits on cr x Br = a{pp ^ q/g —>■ 
X2+^) X Br(x2 7X1) using a Bayesian approach with a flat prior distribution for the signal 
cross section. The resulting upper limit as a function of squark/gluino mass is displayed in 
Fig. ^, for the case where rriq = rrig, along with the theoretical cross sections, as calculated 



using the CTEQ3L parton distribution function (p.d.f.) |jTT|. The hatched band represents 



the range of the theoretical cross sections obtained by varying the supersymmetry parameter 
values with the constraints Br(x2 —>■ IxT) = 100% and m^o — m^o > 20 GeV/c^. The 
intersection of our limit curve with the lower edge of the theory band is at a x Br = 0.38 pb, 
leading to a lower limit for equal mass squarks and gluinos of 311 GeV/c^. 
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FIG. 3. The 95% C.L. upper limit on the it x Br as a function of rrig/g assuming an equal 
squark and gluino mass. The hatched band represents the range of the theoretical cross section 
for different sets of MSSM parameter values consistent with the constraints Br(x2 IXi) = 100% 
m^o > 20 GeV/c2. 



and m,-o 



The effect of light sleptons on the squark and gluino decays is studied by varying the 
slepton mass (m^) in MC. The fraction of events containing at least one X2 increases as the 
mi is decreased. When is varied from 500 GeV/c^ to 80 GeV/c^, the fraction increases 
about 25% for nig = nig = 300 GeV/c^ with equal gluino mass. Sleptons with mass below 
80 GeV/c^ have already been excluded |T^. The increased X2 production increases the mass 
limit by approximately 10 GeV/c^. 

A light stop (ti) will also modify the squark and gluino decays and therefore affect the 
X2 production. We investigate this effect by setting = 80 GeV/c^ which approximately 
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corresponds to the current ti lower mass limit [T^. A 15% reduction in X2 production cross 



section is observed. This reduction lowers the limit for equal mass squarks and gluinos by 
about 6 GeV/c2. 

Following the procedure above, we obtain a low mass limit for gluinos (squarks) to 
be 233 GeV/c^ (219 GeV/c^) when squarks (gluinos) are heavy. Again, these limits vary 
approximately 10 GeV/c^ if ti and/or sleptons are light. 



III. SEARCH FOR DIRAC MONOPOLES 



One of the open questions of particle physics is the existence of Dirac monopoles |T3|,|14 
hypothetical carriers of the magnetic charge proposed by P.M. Dirac to symmetrize Maxwell 
equations and explain the quantization of electric charge. If such magnetic monopoles exist, 
then the elementary magnetic and electric charges {g and e) must be quantized according 
to the following formula: 



g = ^, n = ±l,±2,..., (1) 



where n is an unknown integer. 




FIG. 4. Feynman diagram for 77 production via a virtual monopole loop. 

Dirac monopoles are expected to couple to photons with a coupling constant ag = 
(yf^/47r ^ 34 n"^ which is at least three orders of magnitude larger than the corresponding 
photon coupling to the electric charge (tte = e'^/Air ^ 1/137). Therefore such monopoles 
could give rise to photon-photon rescattering via the box diagram shown in Fig. ^ |jl5|,|16[. 
The contribution of this diagram for pointlike monopoles to diphoton production at hadron 
colliders was recently calculated and shown to be significant even for monopole masses 



comparable to the collider beam energy. 

Since the virtuality (Q^) of most incoming photons in the process of Fig. ^ is small |]1 



the interacting partons scatter at very small angles and therefore escape the detector through 
the beam pipe. Thus, a signature for monopoles at hadron colliders is the production of a pair 
of isolated photons with high transverse energies. This process gives a unique opportunity to 
find evidence for Dirac monopoles or to set limits on the monopole mass. Previous monopole 
searches can be found in Ref. [lT9|,pO|. 
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Despite numerous studies, QED with pointlike monopoles is still not a complete theory. 
For example, it is not clear whether such a theory can be constructed to be renormalizable 
to all orders [|I6|. Also, arguments exist (see, e.g. |^) that Dirac monopoles must occupy 



a spatial volume of radius R ~ 0{g'^/M), where M is the monopole mass, to accommodate 
the self energy implied by the large coupling. In such a theory, it is possible that hard 
interactions of a monopole with photons would be weakened substantially by the effects of 
a monopole form factor. 

The data used in this analysis represent an integrated luminosity of 69.5 ±3.7 pb^^ using 
a trigger which required the presence of an electromagnetic object with transverse energy 
Ej- above 40 GeV. This trigger did not require the presence of an inelastic collision, and 
therefore can be used to select low events typical of the process in Fig. ^. 

The following offline selection criteria are: (i) at least two photons with Et > 40 GeV 
and pseudorapidity \t]'^\ < 1.1; (ii) missing transverse energy in the event I^t< 25 GeV; and 
(iii) no jets with Ej, > 15 GeV and \ri^ < 2.5. The jet veto requirement is used to select 
the low process in Fig. |. The trigger is > 98% efficient for this off-line selection. 

Photons are selected from the identified EM clusters by requiring no tracks pointing 
toward the cluster from any of the event vertices. The overall efficiency for photon identifi- 
cation is (73.0 ± 1.2)% per photon. The total efficiency for diphoton is (52.8 ± 1.4)%. This 
includes the efficiency of the I^t veto (99.0 ± 0.5)% as well as the identification efficiency 
for a pair of photons. The above selection criteria define our base sample which contains 90 
candidate events. 

The main backgrounds to photon scattering through a monopole loop are due to: (i) 
diagrams similar to Fig. ^ with other particles in the loop; (ii) QCD production of dijets 
(jj) and direct photons (j'j) (with jets misidentified as photons due to fragmentation into 
a leading 7r° or rj decaying into a pair of spatially close photons, reconstructed as one EM 
cluster), or direct diphotons (77); and (iii) Drell-Yan dielectron production with electrons 
misidentified as photons due to tracking inefficiency. 



Background (i) is dominated by a virtual ly-loop and has been shown to be negligible [22 



The other two background contributions are estimated from the data. The QCD background 
is determined using the jj event sample collected with a single photon trigger, with the 
jet passing the same fiducial and kinematic cuts as the photon. By applying a jet-faking 
probability, we find the QCD background to be 25 ± 8 events. Direct photon and diphoton 
backgrounds are also included in this estimate. 

The Drell-Yan background is calculated from a sample of dielectron events passing the 
same fiducial and kinematic cuts as the signal sample. Multijet contamination of this sample 
is negligible. The probability for a dielectron pair to be misidentified as a diphoton pair is 
found to be (11 ± 1)% by comparing the number of events in the Z peak in the ee and 
77 samples passing loose kinematic cuts. The Drell-Yan background in the base sample is 
63 ± 7 events. The overall background in the base sample is 88 ± 11 (syst) events, in good 
agreement with the 90 observed candidates. 

To optimize the sensitivity of this search to the monopole loop contribution we apply a 
cut on the scalar sum of the transverse energies of all the photons in the event: St = J2i -E't • 
We vary the St cut threshold (S"™™) in 10 GeV steps to achieve an expected background 
of 0.4 events. Such an optimization is based on the fact that for this expected background 
one has a 67% probability of observing no candidate events in the data in the absence of a 
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signal. In such a case pO[, the hmits on the signal do not depend on the exact background 
value or its uncertainties. The agreement between the observed number of events and the 
predicted background as a function of S"™" is illustrated in Fig. ||. The S*™™ = 250 GeV cut 
corresponds to a background of 0.41 ±0.11 events. We set an upper limit for the production 
cross section of two or more photons with J^Ej' > 250 GeV and \ri'^\ < 1.1: 



a{pp ^> 2-f)\sj,>250 Gov,|,77|<i.i < 83 fb 



(2) 



at the 95% C.L. 




Sj. (GeV) 

FIG. 5. Data and expected background as a function of S™™ cut. Points are data, the up- 
per hatched region corresponds to the QCD background, and the lower shaded region shows the 
Drell-Yan background. The ~ 15% systematic error on the background is not shown. 



Since the data are consistent with the background hypothesis, we can set limits on the 
production of pointlike Dirac monopoles. We calculate the acceptance for the monopole 
signal using a fast MC program that generates diphoton events from a monopole loop ac- 
cording to the calculated differential cross section d^a/dEj^drj'^^dr]'^^ with a subsequent 
parametric simulation of the D0 detector. The overall acceptance for the monopole signal 
is found to be (51 ± 1)%. The acceptance does not depend on the monopole mass for masses 
above the typical photon energy (~ 300 GeV) |1^ . 



The total cross section for heavy monopole production at the Tevatron is given by JT? 



cr{pp ^ 'j'j + X) = 57P 



n 



M [TeV] 



fb, 



(3) 



where P is a spin dependent factor p^j23[| : P = 0.085, 1.39, and 159 for monopole spin of 0, 
1/2, and 1, respectively. The estimated error on this cross section due to choice of p.d.f. and 
to higher order QED effects is 30% |T^. Additional uncertainties are associated with the 



77 ~^ 77 subprocess in Fig. ^ and with unitarity considerations. The coupling constant ag 
is replaced with an effective coupling [|T7| obtained by multiplying ag by a factor {E'' /M^, 
where E"^ is the photon energy, typically 300 GeV at the Tevatron. Both unitarity and 
perturbation theory assumptions are satisfied when this factor is <C 1 1]T3|,|T7 . 

Comparing the lower bound of the theoretical cross section corrected for acceptance with 
the cross section limit set by this measurement, we obtain the following lower limits on the 
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FIG. 6. The curved bands show the low and upper bounds on theoretical cross sections 

for 

monopole spin, S = 0,1/2, and 1. The horizontal line shows the 95% C.L. experimental upper 
limit on the cross section. The arrows indicate the lower 95% C.L. limits on the monopole mass 
at each spin value. 



pointlike Dirac monopole mass (see Fig. 

r 610 GeV for S = 
M/n > I 870 GeV for S = l/2 . 
[ 1580 GeV for 5 = 1 

We note that the effective coupling exceeds 1 and unitarity is violated close to the exper- 
imental bound. For values /M > 1, the cross section will grow more slowly, approaching 
the usual 1/M^ behavior of a QED process [jl8[ which satisfies unitarity. Also, for lower 



monopole masses the effective parameter of the perturbation theory used in the calcula- 
tions becomes too large, and therefore one would expect a non- negligible contribution 



of the higher order diagrams with four, six, etc. photons in the final state. The latter effect 
is, however, largely compensated by our analysis cut on the sum of the photon transverse 
energies; if part of the signal cross section is due to the higher order diagrams, the above 
limits are unaffected. 



IV. SUMMARY 



In summary, we have performed searches for supersymmetry by searching for 
and 7^r+ > 2 jets events and for heavy pointlike Dirac monopoles by searching for pairs 
of photons with high transverse energies. Our data agree with the expected background 
production. Within the framework of the MSSM with the choices of the parameter values 
consistent with Br(x2 7Xi) = 100% and m^o — m^o > 20 GeV/c^, we obtain a 95% 
C.L. lower mass limit of 311 GeV/c^ for equal mass squarks and gluinos. Using theoretical 
calculations fl^ we set 95% C.L. lower limits on the Dirac monopole mass for minimum 
magnetic charge (n = 1) in the range 610 to 1580 GeV, depending on the monopole spin. 
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